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Non-perturbative lattice simulations have shown that there is no electroweak phase transition in 
the Standard Model for the allowed Higgs masses, mn ^ 75 GeV. In the Minimal Supersymmetric 
Standard Model, in contrast, it has been proposed that the transition should exist and even be 
strong enough for baryogenesis up to mn ~ 105 GeV, provided that the lightest stop mass is in the 
range 100. . . 160 GeV. However, this prediction is based on perturbation theory, and suffers from 
a noticeable gauge parameter and renormalization scale dependence. We have performed large- 
scale lattice Monte Carlo simulations of the MSSM electroweak phase transition. Extrapolating the 
results to the infinite volume and continuum limits, we find that the transition is in fact stronger 
than indicated by 2-loop perturbation theory. This guarantees that the perturbative Higgs mass 
bound mn ~ 105 GeV is a conservative one, allows slightly larger stop masses (up to ~ 165 GeV), 
and provides a strong motivation for further studies of MSSM electroweak baryogenesis. 



It is known from studies of primordial nucleosynthesis 
that there is a non- vanishing baryon to photon density ra- 
tio in the Universe, r\ w 1CP 10 (for recent reviews, see [jij). 
It is one of the main challenges of cosmology to under- 
stand how such an asymmetry could come about. Indeed, 
different scenarios for producing r\ > abound. 

Among all the scenarios for baryogenesis, one is unique: 
the last instance in the history of the Universe that a 
baryon asymmetry could have been generated, is the elec- 
troweak phase transition 0] . As such, this is also the sce- 
nario requiring the least assumptions beyond established 
physics. In principle, even the Standard Model con- 
tains the necessary ingredients for baryon number genera- 
tion: anomalous baryon number violation, CP-violation, 
and an electroweak phase transition providing for a non- 
equilibrium environment (for a review, see 0). Once an 
asymmetry has been generated, it must also be preserved, 
and this gives a strict constraint on how strongly of the 
first order the transition must be Q. In fact, the con- 
straint on the strength of the phase transition is the most 
rigorous of the constraints mentioned, since it concerns a 
thermodynamical equilibrium situation after the transi- 
tion, and equilibrium physics is much better understood 
than non-equilibrium physics. 

However, it turns out that on a more quantitative level 
the Standard Model is too restricted for baryogenesis. 
The main reason is that the strength of the electroweak 
phase transition depends on the Higgs mass, and for the 
allowed values m# > 75 GeV, there is no electroweak 
phase transition at all [[|j|] . The existence of the baryon 
asymmetry alone thus requires physics beyond what is 
currently known. 

The simplest extended scenarios that allow for baryon 



asymmetry generation at the electroweak phase transi- 
tion, have a Higgs sector which differs from that in the 
Standard Model. A particularly appealing scenario is the 
electroweak phase transition in the MSSM [6-8]. Indeed, 
it has recently become clear that the electroweak phase 
transition can then be much stronger than in the Stan- 
dard Model, and strong enough for baryogenesis at least 
for Higgs masses up to 80 GeV [9-16]. For the lightest 
stop mass mj lighter than the top mass, one can go even 
up to ~ 100 GeV p7| : in the most recent analysis [Q, 
the allowed window was estimated at to# ~ 75 . . . 105 
GeV, 771^ ~ 100 . . . 160 GeV. In this regime, the transi- 
tion could even proceed in two stages ]i7j |, via an exotic 
intermediate colour breaking minimum. This Higgs and 
stop mass window is interesting from an experimental 
point of view, as well, as the whole range will be covered 
at LEP and the Tevatron @. 

Unfortunately, the statement concerning the strength 
of the electroweak phase transition in this regime is sub- 
ject to large uncertainties. The first indication in this 
direction is that the 2-loop corrections to the Higgs field 
effective potential are large and strengthen the transi- 
tion considerably jL0|| . A further sign is that the gauge 
parameter and, in particular, the renormalization scale 
dependence of the 2-loop potential, which are formally 
of the 3-loop order, are numerically quite significant [fl7| . 
Hence a non-perturbative analysis is needed. 

The purpose of this paper is to study the MSSM elec- 
troweak phase transition with lattice Monte Carlo sim- 
ulations, and to extrapolate the results to the infinite 
volume and continuum limits. Since the MSSM at finite 
temperature is a multiscale system with widely different 
scales from ~ irT to ~ , and since there are chiral 
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fcrmions, the only way to do the simulations in practice 
is to use an effective 3d theory [^9|. This approach con- 
sists of a perturbative dimensional reduction into a 3d 
theory with considerably fewer degrees of freedom than 
in the original theory [21-23], and of lattice simulations 
in the effective theory. The analytical dimensional reduc- 
tion step has been performed for the MSSM in [12-14,17]. 
Lattice simulations in dimensionally reduced 3d theories 
have been previously used to determine the properties of 
the electroweak phase transition in the Standard Model 
in great detail [24-30] . 

In the regime considered, the right-handed stop field 
U plays an important role in addition to the Higgs 
field. The effective 3d Lagrangian describing the elec- 
troweak phase transition in the MSSM is therefore an 
SU(3)xSU(2) gauge theory with two scalar fields 
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Here Df and Df are the SU(2) and SU(3) covariant 
derivatives, and H is the combination of the Higgs dou- 
blets which is "light" at the phase transition point. The 
U(l) subgroup of the Standard Model induces only small 
perturbative contributions |50|, and can be neglected. 

The complexity of the original 4d Lagrangian is hidden 
in Eq. (|IJ) in the expressions of the parameters of the 3d 
theory. A dimensional reduction computation leading to 
actual expressions for these parameters has been made 
in jl7j for a particularly simple case. Let us stress here 
that the reduction is a purely perturbative computation 
and is free of infrared problems. The relative error has 
been estimated in |l^Jl^], and should be < 10%. 

It is prohibitively time-consuming to study the full pa- 
rameter space of Eq. (|j) with Monte Carlo simulations. 
Thus, we only consider a special parameter choice: we 
take a large left-handed squark mass parameter rriQ ~1 
TeV, vanishing squark mixing parameters, and a heavy 
CP-odd Higgs particle (m A > 300 GeV). We fix tan/3 = 
3, corresponding to mn ~ 95 GeV. We then study the 
3d theory in Eq. (|l|), parametrized by the temperature 
T and the right-handed stop mass parameter rhu {rhu 
determines the zero temperature right-handed stop mass 
through toj « (m 2 op — my) 1 / 2 ). The actual expressions 
used for the dimensional reduction are given in j3l[] . 

The philosophy is now that we determine the non- 
perturbative results for the continuum theory in Eq. (Q) 
through lattice simulations, and compare them with 3d 
perturbation theory, employing the same 3d parameters. 
To be more precise, we compare with 2-loop 3d pertur- 
bation theory in the Landau gauge £ = and for the 
MS scale parameter /I = T, values which have been used 
in [fiJI , as well. This allows one to find out whether there 
are any non-perturbative effects in the system. Once this 



has been done, one can go back to a more complicated 
situation and study it perturbatively, adding to the per- 
turbative results the non-perturbative effects found here. 
As the reduction step is purely perturbative, the non- 
perturbative effects found with the 3d approach apply 
also to the effective potential computed in 4d Jl~0| , |T6| , |l8[ | . 

To perform lattice simulations, we discretize the the- 
ory in Eq. (m with standard methods (see 31 1). The 



lattice parameters are expressed in terms of the lattice 
spacing a and the continuum parameters through 2-loop 
relations p2| which become exact in the continuum limit. 

Well controlled infinite volume and continuum limits 
are essential in order to obtain reliable results. Thus, for 
each point in the parameter space, we always perform 
simulations with several lattice volumes and extrapolate 
to the infinite volume. We use the lattice spacings ob- 
tained through 
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= 12,20, 
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9S3 = 9sT ~ T is the 3d SU(3) gauge coupling 
and a is the lattice spacing. The fact that we use just 
two values of 0s, only allows a linear extrapolation to 
the continuum limit 0s — oo. However, it is understood 
analytically that the dominant corrections are linear |33) , 
and moreover, linear extrapolations work extremely well 
for the case of the Standard Model 25 3Pfl- 

All in all, we have performed 42 different Monte Carlo 
runs: combinations of lattice sizes and parameters. The 
total cpu-time was ^7.5 node-years on a Cray T3E. 

The physical quantities we discuss here are the criti- 
cal temperature T c , the scalar field expectation values, 
and the latent heat. Quantities such as the latent heat 
enter, for instance, the estimates for the nucleation and 
reheating temperatures (see, e.g., Q), which are needed 
to decide whether the scalar field expectation values rel- 
evant for cosmology should be taken at T c or some lower 
temperature. 

1. The phase diagram and the critical temperatures. 
The general phase structure of the theory is expected to 
be the following Jl7j . The system has a first order transi- 
tion at T c ~ 100 GeV for m v < 65 GeV. This transition 
is strong even though mjy is large, due to the stop loops. 
As rhu becomes larger {mi smaller), the transition gets 
even stronger, and then at some point one may get a 
two-stage transition. The existence of a two-stage tran- 
sition depends on the parameters of the theory, and for 
large squark mixing parameters the two-stage region is 
not reached Jl8| . 

Our numerical results are shown in Fig. [j]. It is seen 
that the phase diagram is qualitatively the same as in 
perturbation theory, although the critical temperatures 
and the triple point have been displaced by a few GeV. 
We have data at 0s = 20 only at rhu = 50, 65 GeV, 
and the continuum extrapolation is possible only at these 
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FIG. 1. The phase diagram and the critical temperatures. 
The continuous lines are from the 2-loop perturbative effec- 
tive potential in the Landau gauge. Open symbols correspond 
to infinite volume extrapolations, and filled symbols to con- 
tinuum limit extrapolations. 
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FIG. 2. The latent heat. 

points. Nevertheless, we expect similar (small) effects at 
the other points. As of now, we have no clear theoret- 
ical explanation for the discrepancy between the lattice 
results and perturbation theory: the reason might be, 
e.g., a three-loop perturbative effect, or a genuine non- 
perturbative contribution. 

2. Latent heat. 

The main result of this paper is shown in Fig. |^, which 
shows the latent heat. It is the most important gauge- 
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FIG. 3. The scalar field expectation values in the broken 
phases at T c . 



invariant physical characterization of the strength of 
a first order transition. We observe that the non- 
perturbative transition to the standard electroweak min- 
imum at rhu < 67 GeV is significantly (up to 45%) 
stronger than the perturbative transition. In the regime 
™>u <; 67 GeV where there is a two-stage transition, a 
comparison with perturbation theory is more difficult as 
the whole pattern is shifted to the right, but the qualita- 
tive behaviour is the same. 

3. Scalar field expectation values. 

The Higgs field vacuum expectation value Vh is the 
object by which one usually characterises whether the 
phase transition is strong enough for baryogenesis H||, 
the requirement being vh/T>1. As such vh is, how- 
ever, a gauge dependent quantity. If one computes it 
in the Landau gauge (vjj), as is usual, then in terms of 
gauge-invariant operators the same expression would be 
non-local. On the other hand, there is a simple local 
gauge- invariant quantity closely related to vh, namely 
H^H - v 2 H /2. The problem with i?t H is that being a 
composite operator, it is a scale dependent quantity in, 
say, the MS scheme. We hence define on the lattice 
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which is a natural gauge-invariant generalization of 
vfc/T, and can be measured in simulations. Note that 
with respect to 4d units, there is a trivial rescaling by T 
in the H^H appearing in Eq. (0). 

The numerical results for vh /T, vjj /T are shown in 
Fig. [| Again, we observe a value larger than in per- 
turbation theory in the regime % < 67 GeV. More- 
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over, in qualitative accordance with perturbation the- 
ory, there is a rapid increase in vh/T c in the regime 
of the two-stage transition, rhu > 67 GeV. The rela- 
tive non-perturbative strengthening effect is smaller than 
for the latent heat, which is easy to understand since 
L oc A(iftjj) _ Av 2 H [||, implying SL/L ~ 2Sv H /v H . 

In conclusion, at least for the parameter values studied 
(m H ~ 95 GeV, mj ~ 150 . . . 160 GeV), the electroweak 
phase transition is significantly stronger than indicated 
by 2-loop perturbation theory. This implies that the pre- 
vious perturbative Higgs and stop mass bounds for elec- 
troweak baryogenesis are conservative estimates. In par- 
ticular, the electroweak phase transition could be strong 
enough for baryogenesis for all allowed Higgs masses 
in this regime (m# < 105 GeV) [Q. Due to the non- 
perturbative strengthening effect seen, the stop mass 
could be slightly larger than the perturbative value, up 
to mi ~ 165 GeV. For the smallest stop masses, on the 
other hand, there is the possibility of a two-stage tran- 
sition, in which the Higgs field gets an extremely large 
vacuum expectation value. 

These results provide a strong motivation for precise 
studies of the non-equilibrium CP-violating real time dy- 
namics and baryon number generation at the MSSM elec- 
troweak phase transition. 
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